CERN-PH-TH /2008-023 
IFJPAN-IV-2007-12 
6 March 2008 



Gauge invariant sub-structures of tree-level 
double-emission exact QCD spin amplitudes 

A. van Hameren0 

IF J- PAN, ul. Radzikowskiego 152, Krakow, Poland 

and 
Z. Wa^ 

IF J- PAN, ul. Radzikowskiego 152, Krakow, Poland 

and 

CERN, 1211 Geneva 23, Switzerland 



Abstract 

In this note we discuss possible separations of exact, massive, tree-level 
spin amplitudes into gauge invariant parts. We concentrate our attention on 
processes involving two quarks entering a color-neutral current and, thanks 
to the QCD interactions, two extra external gluons. We will search for forms 
compatible with parton shower languages, without applying approximations 
or restrictions on phase space regions. Special emphasis will be put on the 
isolation of parts necessary for the construction of evolution kernels for in- 
dividual splittings and to some degree for the running coupling constant as 
well. Our aim is to better understand the environment necessary to opti- 
mally match hard matrix elements with partons shower algorithms. To avoid 
complications and ambiguities related to regularization schemes, we ignore, 
at this point, virtual corrections. Our representation is quite universal: any 
color-neutral current can be used, in particular our approach is not restricted 
to vector currents only. 
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1 Introduction 



One of the essential steps in the construction of any algorithm for multi-particle final 
states is the appropriate analysis of the phase space parametrization. In the PHOTOS 
Monte Carlo [1] for multi-photon production, such an exact phase space parametrization 
is embodied in an iterative algorithm, the details of which are best described in [2], but 
the control of the relative size of sub-samples of distinct numbers of final state particles 
requires a precise knowledge of the matrix elements including virtual corrections. In the 
KKMC Monte Carlo, the approach to phase space generation is different, but the necessity 
to control matrix elements and their truncation is again essential [3,4]. 

Iterative procedures for parts of amplitudes, which are at the foundation of exponentia- 
tion [4,5] and structure functions [6-10] were exploited for the sake of use in KKMC Monte 
Carlo. In particular the description of the dominantly s-channel process e^e~ z/gPeTT 
where i-channel VT-exchange diagrams with gauge boson couplings contribute to matrix 
elements provided an interesting example [11]. These studies were motivated by practi- 
cal reasons, but also pointed at quite astonishing properties of trcc-lcvcl spin amplitudes, 
namely that they can be separated into gauge invariant parts in a semi-automated way, 
easy to apply in the Kleiss- Stirling methods [12, 13]. 

One could ask the question whether similar techniques can be used in QCD, whether 
they are of any practical use, and in fact to which degree they were already included in 
previous work, such as for example at the foundation of parton shower algorithms (new 
or well established) [14-22] or for other, fundamental or phenomenological applications 
[23-30] . The common idea between all of these papers is to separate approximate or exact 
results into parts, often with the help of iteration. 

We will not elaborate on these points here, this broad subject by far exceeds the purpose 
of the present paper. We only want to mention the existence of possible limitations in such 
strategies, if applied to parton shower applications [31]. Our study will be purely technical, 
and we will also ignore any possible complications related to confinement for example. As 
was done for QED, we will study possible ways of separating amplitudes, but now for the 
process of two quarks and two gluons entering a color-neutral current. First, we try to 
identify sub-structures of the amplitude which would be useful for possible applications. 
Then we would like to verify to which degree they can be used to localize parts of the 
amplitude related to eg. evolution kernels and the running coupling constant. We expect 
such expressions to be identifiable at least in approximations valid in certain regions of 
phase space, dominant for specific purposes, like in applications using pt ordered phase 
space, but we hope to localize them partly already at the level of exact amplitudes. 

However, we also expect that the choice of what we call subtraction terms will be more 
ambiguous than in the QED case. It concerns the terms which are subtracted at one place of 
an expression, and added at another one. For QED, these terms could be constructed from 
parts of lower-order amplitudes, and resulted in a semi-automated method of organizing 
amplitudes. There, the most singular parts of eg. double photon bremsstrahlung amplitudes 
could be deduced from amplitudes for single photon emission. In fact, those properties are 
known since long and necessary for exponentiation. In the calculation of the so-called /32 
for double-photon emission, the most singular parts are subtracted from the amplitude. 
This is not only convenient but also necessary, since those terms are already included, in a 
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rigorous way, through the exponentiation of lower order amphtudes. In our present study, 
it will not be possible to construct the subtraction terms in such an unambiguous way. Our 
expressions concern the amplitude for the process qq Jgg, but because they are exact, 
the incoming and outgoing particles can be easily interchanged to describe for example the 
process qg Jgq. The amplitudes for the two processes are basically the same, but the 
lower order processes, for which they may be considered as real-emission corrections, are 
different. This clearly points at ambiguities if the form of lower order amplitudes is to be 
used for the definition of subtraction terms. Obviously, by inspection of the equations for 
the evolution of the parton showers or PDFs, both of the two processes are potentially of 
similar phenomenological importance. 

Our paper is organized as follows. In Section [2] we present our notation and general 
strategy to organize the amplitude. The treatment of different possibilities for this or- 
ganization concerning the color part of the amplitude is distributed over Section [3] and 
Section m We start with a formulation close to the one used in previous papers on QED, 
moving consecutively towards shorter and more natural forms for QCD applications. In 
Section [5] we explore certain limits, in which some parts of the amplitude can be dropped 
out. The summary in Section [6] closes the paper. 

2 Notation and general strategy 

The exact spin amplitude for the process qq —>■ Jgg can be written as an expansion in com- 
binations of the SU(Ai"c)-generators, for example the combinations {T",T^} and [T",T'']. 
Here, T" is the color generator associated with gluon number 1, which has momentum ki 
and polarization vector ei, and is the color generator associated with gluon number 2, 
which has momentum k2 and polarization vector 62- Another option would be to use, for 
example, the combinations T°'T'' and T^T"". But let us stick to the former example for 
now, and return to this issue later. The amplitude can then be expressed as 

M"^^ = ^ v{p) {[T\ T^ti'^l + {T\ T^}/{1'2} j u{q) , (1) 

where v{p) and u{q) are the spinors associated with the anti-quark, which has momentum 
p, and the quark, which has momentum g, respectively. We do not specify the spin or 
color states for the quark fields, any choice can be used. This type of separation of the 
spin amplitude into gauge invariant parts is known [32] and exploited since long time. 
The main task is now to calculate the coefficients 7^^'^' and J^^'^^. They can be written 
as sums of terms, each of them containing factors consisting of scalar products of the 
aforementioned four-vectors and contractions of these four-vectors with the Dirac gamma- 
matrices. Furthermore, each term will carry a factor /, associated with the color-neutral 
current J. We want to stress that all expressions we present are valid under the condition 
that the current / is constructed from color neutral objects like {v + 07^), 7^^, ^1, etc., 
although our notation might seem to indicate stronger limitations. 

The first step in this calculation is the construction of the relevant graphs, depicted in 
Fig.[Tl from the Feynman rules. Next, we choose to replace explicit mass terms by fermion 
momenta using the Dirac equation. It is not difficult to extract expressions of the type 
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Figure 1: Feynman graphs for the process qq — > Jgg. 



mentioned before for /l^'^l and /^^'^^ at this point. These expression will contain terms 
with products of zero, two or four gamma-matrix contractions (or "slashed" four- vectors) , 
excluding the object /. Thanks to properties of the gamma- matrices, we can now choose 
how these products are ordered, at the expense of the appearance or disappearance of 
certain scalar products of four-momenta. We choose to order them such that factors /^i,2, 
^1,2/ or ^1,2^2,1, ^1,2^2,1 are not present. This is for the convenience of possible use of the 
Kleiss- Stirling method for coding final expressions similarly as in [3] . 

The game is now to compactify expressions, or organize them in an easily interpretable 
way, by factorizing certain sums of terms. The most important guideline we use is that each 
term, consisting of factors which consist themselves of sums of terms, should be manifestly 
gauge invariant. In this process, we allow for the introduction of what we call subtraction 
terms, that is terms subtracted at one place of the complete expression, and added at 
another place. They can help to make certain factorizations possible, which look desirable 
but for which some terms appear to be missing. 

Such subtraction terms were also used in our earher QED-studies of the same type [11]. 
In fact, the terms were constructed from parts of lower-order amplitudes and were an 
important element in the attempt to define a semi-automated method to obtain expressions 
for spin amplitudes also used in Monte Carlo programs. For the two-photon case, the 
following subtraction terms were used: 



and similar terms with q replaced by p. As mentioned before, these terms are constructed 
from parts of the one-photon amplitude, which is given by 






(3) 



M'' = v{p)T''I^^^u{q) , 



(4) 



with 
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Note that each segment encapsulated by square brackets is manifestly gauge invariant. 

As in the case of QED, we will use subtraction terms ([2]) and ([3]). In addition however, 
we will also use 

^{1,2} _ _ 1 J h-h ^1-62 fc2-ei 

and a similar term with q replaced by p as subtraction terms. Note that, in contrast to 
([2]) and ([2]), subtraction term ^ is not constructed from elements present in the single 
bremsstrahlung amplitude given in Eq. ([5]). In particular, it introduces an artificial singu- 
larity l/(/ci-/c2) which would be worrisome for QED applications. We, however, are not 
concerned with that, since it points at the singularity of an intermediate gluon virtuality, 
with which we have to deal anyway in a QCD amplitude. 

We already mentioned the examples of an expansion in \T'^,T^] and {T",T^}, and an 
expansion in T"T^ and T^T"". We will deal with both cases in the next sections. We 
start with the first alternative because the {T°, T^}-part has exactly the same form as a 
two-photon amplitude in QED. The second part, proportional to [T'*,T^], is new in that 
respect. Naively, one would expect it can, to a large degree, be identified as a correction to 
the single gluon emission - that is the case when this emission becomes virtual and decays 
into a pair of gluons. For that part of the amplitude we will use subtraction terms quite 
similar to the ones discussed here. 

The second alternative, of an expansion for the color contents in terms of T"T^ and 
T''T"', is inspired by approximations exploited for example for parton showers and in par- 
ticular dipole models. Then, the gluons can be considered to appear consecutively (in one 
sense or another). 

The most compact expressions for the amplitude will be obtained with a mixed choice, 
in which all four combinations [T",T^], {T'^,T^}, T'^T^ and T^T'^ of color generators are 
used. The coefficients for each term is not unique in this case of course, but we will be 
able to find a choice which seems to be optimal. 



3 QED-inspired picture 

In this section, we investigate expressions for the two-gluon amplitude obtained when 
expanded in its color contents as in Eq. ([T]). We start our analysis by collecting terms 
proportional to the anti-commutator of the color generators. It can be considered the 
less complicated one, because it does not include terms originating from the triple-gluon 
vertex. Note that it is identical with the expression for QED amplitudes of the process 
e~^e~ — > i^fMi^^j'y as described in paper [11]. This is possible, because the anti-commutator 
projects out the triple-gauge coupling. The choice of the subtraction terms affects the final 
form of the results. In case of QED we could limit ourselves to the gauge invariant parts 
available from the single photon (gluon) amplitude, leading to a seemingly unique result. 
The case of QCD is more complex, and the choice of subtraction terms is rather moti- 
vated by aestetical considerations, for example the identification of simple gauge invariant 
expressions which depend simultaneously on the momenta p and q, leaving more lengthy 
expressions dependent on only one fermion momentum for further manipulations. In fact, 
for the anti-commutator part, we recall the QED result, continuing the separation into 



5 



parts somewhat further, hoping for new parts with properties which were of no interest 
for [11] or even inconsistent with the structure of QED singularities. 



3.1 {T^T^j-part 

The part of the amphtude proportional to {T",T^} can be represented as sum of twelve 
individually gauge invariant parts: 



J{1,2} 



where 
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(7) 



(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 



^q-ki ki-k2jq-k2 \q-k2 ki-k2jq-ki 

This form should be compared with the one of formula (27) in [11], where the amplitude 
is separated into seven terms A^i to M.^. Our expression ^ equals M.t. Expression 
([9]) added to ffTOj) equals A^g- Expression f|T2l) and f|T3|) equal M.i and M.2 respectively. 
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Expression (fTTj) . equals M.^. The sum of expressions (JHl), (fT5|) and (fTSj) equals A^s. 
Analogously, the sum of expressions (fT6l) . f|T7|) and (fT9l) equals A^4. 

Subtraction term ([2]) is integrated with a positive sign into ([H]), and with negative sign 
into expressions fITBl) . (11 7p and (1191) . The analog of ([2]) for momentum p is integrated into 
expression ([8]) with a positive sign, and into expressions f|T4|) . f|T5l) . and f|T8|) with a negative 
sign. 

Subtraction term ([3]) is integrated with a positive sign into iQ and ffTOj) . and with a 
negative sign into expressions flTB]) . f[T7|) and flTI?]) . The analog of ([3]) for momentum p is 
integrated with a positive sign into and ffTU]) . and with a negative sign into expressions 
([HD, ([ISD, and ([IHD. 

Note that contrary to [11] we have allowed here a new, from the point of view of 
fermionic spin amplitudes and Kleiss-Stirling methods [12] ugly, object 61-62. It will appear 
anyway in the commutator part discussed in the next section as a consequence of the triple 
gauge coupling. We have found it in (ITBl) and (IT7I) because of the subtraction term (jS]) 
incorporated with the opposite sign to f|T4l) and f|T6|) . 

As mentioned above, we separated the expressions and 7Vl4 of [11] into several 
new parts, partly because of the new subtraction term. Such a separation is inappropriate 
for QED applications, since the new gauge invariant parts, even though compact, are 
more singular in the soft-photon limit than their sum, and would form an obstacle to 
exponentiation. For QCD, however, this is not so much an issue because of the triple- 
gluon vertex; such a separation will be useful later. That is why we allow the structure 
of apparent singularities to be larger than necessary. As a consequence, the hierarchy of 
dominant and sub-dominant terms is richer, and potential interest in separation into more 
fine parts may go further. 

In [11], the grouping of gauge invariant segments was prepared such that the terms for 

(3^ and of exclusive exponentiation could be identified. Such an approach can be 
productive only in combination with the discussion of a well defined calculation scheme 
like QED exponentiation, or at least an algorithm for phase space generation, in particular 
its pre-samplers and phase space Jacobians. Therefore, the discussion of such an optimal 
choice of grouping would be far too pre-mature here, and we limit ourselves to techni- 
cal properties of the amplitude. Only when similar work for other multi-leg processes is 
completed as well we may return to this issue. 

3.2 [T^ r^]-part 

For the commutator part, we introduce subtraction terms mainly in order to separate 
terms depending simultaneously on both momenta p and q into simple and intuitively 
clear expressions. The part of the amplitude proportional to [T°,T^] is again represented 
as a sum of individually gauge invariant parts and reads: 

r[l,2] _ r[l,2] , r[l,2] .[1,2] .[1,2] .[1,2] .[1,2] .[1,2] .[1,2] .[1,2] .[1,2] 

—-'1 "'"-'2/ "T -'2r "T -'3 "'"-'4/ "T -'4r -'5M -'siB -'SrA -'SrB 

.[1,2] .[1,2] .[1,2] .[1,2] .[1,2] .[1,2] , . 

^ ^61 ^ ^ ^riA ^ ^rW ^ ^7rA ^ ^7rB K'^^ J 
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where 
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1^1 thh _^ ^2|^2 fhii 
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This time, we do not have a good motivation for the particular form. We mainly tried to 
keep it analogous to I^^''^^ term by term. Only the terms with a subscript starting with 7 
do not have an analogue in The subtraction terms 

Ki-K2 \P K2 Q K2 / 

4f ^-^(Mi^i^Mii:^) (39) 

Kq-h q-h q-ki q-k2) ^ ' 

o[i,2] ^ 1 n 1 ( 2g- 62^2-61 2g-eifcre2 \ 

''^ 2U-ki + q-k2-kr-k2\ q-k2 q-ki J ^ ^ 

are quite analogous to the ones used in the previous subsection and we will not elaborate 
on them further. The S^^^p \ 82^^^ and S^^^^^ can be obtained analogously. 



4 Picture of consecutive gluon emission 

Encouraged by the rich substructure of gauge invariant parts we found in our expressions, 
we search for other possible organizations of terms. In particular, we try to incorporate an 
ordering of the gluons, which has proved its usefulness in a broad spectrum of applications 
and calculations concerning QCD. 



4.1 Color-ordered amplitudes 

Instead of an expansion in the commutator and the anti-commutator of color generators, 
which is useful for the comparison with existing QED results, let us simply use generator 
products to express 

The so-called color-ordered amplitudes are obtained from the coefficients /(^'^) and /(^'^) 
by including the spinors v(p) and u{q). Thanks to this ordering the expressions shorten 
significantly. For the T'"T'''-part, we find 
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(46) 
(47) 
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The part proportional to T^T" is obtained by a permutation of the momenta and polar- 
ization vectors of the gluons, and reads 

\p-k2 ki-k2 2p-A;2/ \'2q-ki k2-ki q-kij 

^ P-ki f _ ki-e2 _ kh \ f 1^ _ ^2-61 _ kh \ f /^gx 
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j^jfl \ / k2-ei ki-e2 62-61 \ 
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(53) 



4 q-k2 + q-ki - k2-ki\ k2-ki 
Each line of the above expressions is individually gauge invariant. 

4.2 Mixed representation 

As we will see now, the expression for the amplitude obtained with a mixed, over-defined, 
basis will not only be shorter than the previously discussed ones, but will also lend itself 
better to analyze the properties of the amplitude. Thanks to such a representation not 
only the exact results, but also expressions dominant in some regions of the phase space 
(obtained simply by truncation) will be more compact. The decomposition of the amplitude 
looks as follows 

= \v{p) (t^T'I^^^ + T'T'^lti! + r , + {T^ T'}I^2i^)u{q) , (54) 

with the coefficients given by 

^(1,2) _ ^ P^ _ ^2-61 _ kh ^ ^1-62 ?-62 



p-ki k2-ki 2p-kiJ \2q-k2 ki-k2 q-k2 

_^ Pv^2 f p-ei k2-ei kh \f p-^2 ^1-62 ^2^2 1^ /^^x 

p-ki+ p-k2 - ki-k2 \p-ki k2-ki 2p-ki)\p-k2 ki-k2 2p-k2/ 



_^ f q-ki f q-ei /c2-6i ^iii\(q-e2 /cr62 ^25^2 



q-ki + q-k2 — ki-k2 \q-ki k2-ki 2q-kiJ\q-k2 ki-k2 2q-h 



j{2,i) _ ^ p-e2 _ ki-e2 _ kh ^ hh ^2-6i q-ei 



p-k2 ki-k2 2p-k2 J \2q-ki k2-ki q-ki 

_^ Pvfcl / p-62 fel-62 ^2|^2 \/ P-61 ^2-61 fl^l . 

p-k2 + p-ki — k2-ki \p-k2 ki-k2 2p-k2) \p-ki k2-ki 2p-ki/ 

_^ J y-/,-2 f q-e2 ki-e2 hi'2\fq-ei k2-ei I^i^i 



q-k2 + q-ki — k2-ki \q-k2 ki-k2 2q-k2j\q-ki k2-ki 2q-ki 
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-'mix — ^ / 1 , . , , , + 



2 \p-ki + p-k2 - ki-k2 q-ki + q-k2 - ki-k2 J \ki-k2 ki-k2 ki-k2 
1 1 [hhhh-hhhh 



4 p-ki + p-k2 - ki-k2\ ki-k2 
1 1 / hhhh - hhhii 



(57) 



4 q-ki + q-k2 - ki-k2\ ki-k2 

r{l,2} 1 1 , 1 \ ^1-62 ^2-61 



^'"'^ 2\p-k^+p-k2~k^-k2'^ q-k^ + q-k2-kvk2)\ k,-k2 ^i-^^J . (58) 

The choice for the above representation is not unique, and at this point it seems to be 
based mainly on aesthetic grounds. Again each hne of the above expressions is individually 
gauge invariant. In order to justify our choice a bit further, we change our strategy in the 
construction of the expressions. So far, the only manipulations we performed consisted of 
the reorganization of terms, and the introduction of our so-called subtraction terms. All 
these terms were written explicitly in terms of momenta and polarization vectors. Now, we 
will leave this path, and start to introduce new objects useful to compactify the expression 
even further. The objects are the momentum of the virtual gluon 

f^i+2 = ki -\- k2 , (59) 

and the four- vectors 

^ _ fc^-fc£/ei^W^_ A _ k^-k^ {Tiji^ihhh) .... 

2k,-k2\ krk2 '''■''V ' ''^+^"2fcrfc2 k,-k2 ' ^ ' 

which represent the effective polarizations of the virtual gluon piled together with its 
propagatorjfl Notice, that the expressions are gauge- invariant. We immediately see that 
the first line of Eq. (IFr|) is equal to 



p-ei+2 g-ei+2 



^p-ki+2- ki-k2 q-ki+2-ki-k2 
Furthermore, we will prove in the next section that 



(61) 



2kyk2 

so that we can write 



l+2 + i77l+2)^l+2 , (62) 



.[1,2] ^ ^{ P-ei+2 g-ei+2 

mix 



p-ki+2 - ki-k2 q-ki+2 - ki-k2 
1 1 



2 p-kij^2 -ki-k2 



i+2 + i7Vi+2)ft+2/ (63) 



•^Proper units of energy will appear only after including the appropriate factor from the phase space 
Jacobian (this is usually done, only after some assumption is made on phase space regions to be considered). 
That is why wrong units for this polarization vector are not as disastrous as they may seem. Note also, 
that ei+2 is parallel to ei+2 and has a significant component along the fci+2-direction. This is another 
reason to be cautious with their physical interpretation. 



11 



Notice the similarity of this expression with Eq. ([5]) for the single-gluon emission. In the 
next section, we will show that in the limit when ki-k2 becomes zero, whether it be a soft 
or a coUinear limit, the quantities 

ei-/c2e2-fci ^ ^ _j ^ Tr(7^^i|^i^2|^2) 



61-62 and ; — ; (64) 

fci ■ fc2 / ki- k2 

stay both finite. In particular, when the polarization vectors 61 and 62 become parallel, 
61+2, absent anyway in the first line of ([M]), becomes zero, indicating that it contributes 
negligibly compared to 61+2, and making the fl63|) even more similar to ([5]). Also, 

P Q \ ki — ko , ^ 



p-ki + p-k2 - ki-k2 q-ki + q-k2 - ki-k2 J ki-k2 

remains finite when ki and k2 become collinear. This property of the expression present 
in the first line of ( 157|) is a technical manifestation of the absence of longitudinal gluons, 
but this time for the virtual gluon. Again, the similarity to ([5]) manifests itself. There, 
the cancellation is exact, and is a consequence of gauge invariance in combination with the 
fact that the gluon is on mass-shell. 

As mentioned, the proof of Eq. (1621) will be given in the next section. Here, we just 
want to mention that 

= 4^r^^ + 8k,. k2 k,.k2 • ^''^ 

4.3 Towards an interpretation at the exact amplitude level 

The expression for the amplitude can be compactified further with the help of the four- 
vectors ^ ^ 

«i = ei - fci and = 6^ - ^^—^ k^ . (67) 

ki ■k2 ki- k2 

These manifestly gauge-invariant objects satisfy 

ai-ki = ai-k2 = a2-ki = a2-k2 = , (68) 

and can be understood as a particular choice of gauge for 61, 62. Alternatively the terms 
ki ei-k2/ki-k2 and fcg 62-^1/^1 -^2 can be understood as defining subtraction terms. The 
coefficients of the amplitude flM|) can now be written as follows: 



7-(l,2) 
mix 



p-ai (^1^1 \ / ^2^2 q-a2 



p-ki 2p-kiJ \2q-k2 q-k- 



^2 



_^ / p-ai jxji \ ( p-a2 j ^ (^gg^) 

p-ki + p-k2 - ki-k2 \p-ki 2p-kiJ\p-k2 2p-k2; 



q-ki f q-ai \ f q-a2 |^2(^: 



q-ki + q-k2 — ki-k2 \q-ki 2q-kiJ\q-k2 2q-h 



*2 
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j(2,i) _ ( p-a2 ^2f2 ^ /f( fifi q-ai 



p-k2 2p-k2j \2q-ki q-ki 

_^ f _ \ f P:^ _ y^ifi j ^ ^^^^-^ 

p-k2 + p-ki - k2-ki \p-k2 2p-k2j\p-ki 2p-ki ' 

^ q-k2 [q-a2 hh\ [ Q'^'i hh 



q-k2 + q-ki - k2-ki \q-k2 2q-k2j\q-ki 2q-ki 



r[i.2] _ ^ /rf P-ki -p-k2 q-k2-q-ki ^ f ai-a2 

-'mix — ^ , . , , , + 



2 \p-ki + p-k2 - ki-k2 q-ki + q-k2 - ki-k2 J \ ki-k2 

11/ ^ll^l42J^2 — <li'2f-2ilf-l 



4: p-ki + p-k2 — ki-k2\ ki-h 
1 . 1 



^2 



(71) 



4 q-ki + q-k2 — ki-k2\ ki-k 



^2 

= ~l Kp-h+pl2~k,.k2 + g-fci + g-L-fci-fcJ ^ ~ ""'""'^ ^^^^ 

The four- vectors ai, a2 appear to be particularly useful for the proofs promised in the 
previous section. Let us note first that in terms of these four-vectors the "polarization 
vectors" of the virtual gluon defined in f l60ll are given by 

k^-kl fcf-fc,^ iTr(7Vi^i^2^2) 

and that 

To proof Eq. (!62|) . we will first show that the formula holds in a particular Lorentz- frame. 
Since the expressions are manifestly Lorentz invariant, that is sufficient for the formula to 
hold in general well. Let us choose a rest-frame of the gluon-pair in which ki and 

^2 are back-to-back and ai points along the a;-axis: 

fcf = A;(l, 0,0,1) , A:2^ = A;(1, 0,0,-1) , < = (0,1, 0,0) , = (0, c, s, 0) , (75) 

with 

k = \/h^h/2 and + = 1 . (76) 

In our frame, we have 

iiM2h - (f^2M4i = -{hh - Mi)c - 4:k'^i^-f^s . (77) 

Observing that '~f^'~f'^ = 17^7%'^, + |i^2 = 2/^7° and that ^1 — ^2 = —2/^7^, we rewrite the 
above expression as 

^2^i)(c + i7's) . (78) 



If 1^2f2 - ^2f2^1f 1 = -If lf2 



We use the definition 7^ = i7°7^7^7'^. 
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On the other hand, we have 



{'^h-k2)h+2h+2 = ihh - hh){-^i-^2) = ~{hh - hh)c , (79) 

and ^ 

{2h-hWi+2h+2 = ihh - hh f^'^^/f'^'^'^ = -ihh - hh> ■ (80) 

Now, it is straightforward to complete the proof of ( l62l) by ehmination of the (frame 
dependent) expressions on the r.h.s. of the above equations. 

For the investigation of the behavior of flM|) in the hmit in which ki and /c2 become 
colhnear, we should not use a Lorentz frame in which these momenta are forced to be 
back-to-back. We keep the three- vector of ki along the 2;-axis but give the three-vector of 
/c2 a component along the y-axis: 

A;r = A;?(l, 0,0,1) , k^ = k',{l,0,y, z) (81) 

with 

y-' = l-z^ = {l-z)il + z) . (82) 



The collinear limit corresponds to the limit z ^ 1. According to (l68i) . the four- vectors ai 
and a2 must have the form 



with the only further restriction that 

c\ + sl = -ei-ei , C2 + S2 = -62-62 . (84) 



Explicit calculation gives 

'61-^2 62-fcl 



ki-k2 

and 



- - 61-62 ) = -ai-a2 = C1C2 + S1S2 , (85) 



iTr(7Vi^i^2^2) \TT{^'iM2h) _ 



ki ■k2 ki- k2 

and we see that both quantities are independent of z 



4.4 Properties of the amplitudes 

Before we discuss the properties of the expressions under certain kinematical constraints, 
eg. orderings, let us analyze here the building blocks which appeared already now, at the 
exact amplitude level. For that purpose we will exploit the form of the amplitude as given in 
Section |4?T] and Section W?2\ In principle, such organization remains valid independently of 
whether there are gluons or fermions in the initial state. Nonetheless, we will progressively 
start to focus our attention on the case where the incoming states are fermions and carry 
dominant momenta; p and q respectively. 
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In Section 13.1 of [33] it was shown in a pedagogical manner, that for example in case 
ki becomes collinear with p the factor 



gives, after phase space integration over the appropriate region, the Altarelli-Parisi kernel. 
It can thus be understood as its precursor at the spin amplitude level. Like this, of course, 
such an expression makes no sense, since it is gauge-dependent. Only a gauge invariant- 
object like 



can be interpreted unambiguously. It can be understood as a precursor at the spin ampli- 
tude level for the emission of ki from a dipole-source spanned by p and k2, but this time 
including the effect due to the fermion spin as well. Not only this assures gauge invariance, 
but also gives a scale for the appropriate logarithm if the integration would be completed. 

At this point, it seems straightforward to interpret expressions like fl88p present in (ISSj) 
or (1561) as factors for the emission from incoming fermions and with angle distributions 
controlled by the direction of the other gluon, completing at the same time the second 
arm of the dipole. This interpretation is based mainly on the visual appearance of the 
expressions, ignoring the structure of singular terms and/or interferences, and therefore 
might be misleading. On the other hand, we want to note that our observations are similar 
to the ones which can be found in the literature, where indeed such interpretation was 
performed for the sake of constructing parton shower models, even though usually some 
approximations were used then. 

Let us look at (l88l) . as present eg. in the first line of (!55|) . Independently whether we 
interpret this factor as the description of the emission of a gluon or not, the momentum 
entering J from the left side is p — ki. If the collinear limit can be used, and thus to a 
good approximation J varies slowly with an eventual virtuality {p — ^i)^, then the effective 
intermediate state of a fermion of momentum p — ki can be used to simplify the description. 
A similar line of arguments is used in [33] (and since long in any formulation of the PDF 
evolution) and we will not repeat it here. Note that the intermediate effective gluon 
definition introduced in Section is valid all over phase space, even in phase space regions 
where it does not have any physical meaning anymore. We do not have such descriptions 
for an intermediate effective fermion available at this moment: we would have to stick to 
the low virtuality limit, and thus have to introduce an approximation. 

Let us shift our attention to the second line of (155|) . For the simplified picture of the 
previous paragraph to work, one of the effective incoming momenta must be a collinear 
projection of p — ki — k2 (or q ~ ki — ^2), while the other one simply remains q or p. No 
ambiguity appears in this respect. Regarding the factors {q-ki) / {q-ki + q-k2 — k2-ki) and 
{p-ki) / {p-ki +p-k2 — k2-ki) we want to mention that, in case of collinear configurations, 
they are indispensable for the redefinition of the spinor normalization used in the definition 
of evolution kernels. As we do not want to limit ourselves to such configurations, but want 
to use expressions valid all over the phase space, definitions of spinors for intermediate, 
seemingly on mass-shell, fermion states should remain valid everywhere as well. We do not 
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have such definition available. We want to conclude the discussion about the interpretation 
of the emission-like factors with the remark that outside of some specific regions of phase 
space, the interferences between contributions of different chains (as given by eg. fl48l) ) can 
be of the same order as their squares, and thus cannot be ignored. 

Let us now turn our attention to fl571) . It is rather tempting to interpret its p- and q- 
dependent parts as the real-emission contributions to the running of the coupling constant 
(of the single gluon emission). The most singular parts of these contributions cancel each 
other partly as we can see in Eq. fl^S]) . Such a partial result may already provide a hint 
on the possible optimal choice of scale to be used as an argument of the running coupling 
constant. 

The role of the remaining parts of the amplitude, expression (!58l) . is less clear, although 
it resembles the running coupling constant part. It may equally well be interpreted as a 
genuine second-order part of the amplitude, which can not be interpreted in the language 
suitable for resummation at all. Note that, from the point of view of pq kik2J kinemat- 
ics, this contribution is less singular than any of the previously discussed ones. 

In the next section, we verify how our somewhat naive interpretation will fare in regions 
of phase space where iterative descriptions such as BFKL [34,35], DGLAP [6,7,36] or 
CCFM [37-39] are valid and are known to be useful. 

5 Picture of ordered gluon emission in dipole lan- 
guage 

In the previous chapters we have collected several forms for exact spin amplitudes. In 
the following subsections, we concentrate on cases of special limits, usually associated with 
some type of ordering. We will attempt to define amplitudes which are dominant in certain 
regions of phase space, but nonetheless valid all over phase space, and which differ from 
the complete amplitude only by gauge invariant and analytically available expressions. In 
every presentation of the following subsections, we will start from the exact spin amplitudes 
as given in Section 14.11 and 14. 2[ but of course our aim is to define correspondence between 
our exact fixed order result and approximate ones, which are at the heart of successful 
Monte Carlo programs and are valid at higher orders as well. 



5.1 X ordering and soft gluon limit (BFKL) 

Let us restrict our attention to the region of phase space where y/s ^ k2 and /c^A;^ — 
ki-k2 (all in the reaction rest frame). As a consequence, we can use the conditions 
p-ki 3> p-fc2 S> ki-k2 and/or q-ki ^ q-k2^ ki-k2 for the localization of dominant terms. 
Such a choice is consistent with the BFKL approximation [34,35]. Under these constraints, 
the T"T^-part flS^ of our expression reduces to 

.(1,2) ^ J P^ei _ kr_ei \ f fcrCa _ J q^ _ k2-ei \ / _ ki^\ 

\p-ki k2-ki)\ki-k2 q-k2) \q-ki k2-ki)\q-k2 ki-k2j ' ^ ' 
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All other parts are smaller. For the T T°-part (156|1 we get 

.(2,1) ^ J _ kreA f k2-ei _ / _ ki^\ f _ ^fc^ei \ 

'"^ \p-k2 ki-k2j\k2-ki q-kj^\p-k2 k,-k2j\p-k, k2-kj ' ^ ' 

After some short manipulations, and since the contributions from Eq. flFTI) and Eq. fISS]) are 
also negligible thanks to the conditions, the full amplitude reduces to: 



J^fpRL = \ v{p)fu{q) 



rparpb 



q-ei p-ei\fq-e2 ki-e2 



q-ki p-kij\q-k2 ki-k 



^2 



p-ei q-ei\fp-e2 ki-e2 



p-ki q-kij\p-k2 ki-k 



(91) 



A picture of linked dipoles is manifest in this expression. It is also worth mentioning that 
( |57|) . suspected to contribute to the running of the coupling constant, can be neglected 
under our conditions. The obtained amplitude (19T!) . is a well defined part of the exact one, 
and at the same time is consistent with the BFKL approximation. Note the absence of 

fci-e2 k2-ei 

'^k;^2h^2 ^ ' 

in (1911) . It is a clear manifestation of large destructive interferences between gauge in- 
variant parts of our expressions for the amplitude given in Sections 14.11 and 14.21 A crude 
distribution, of some hypothetical Monte Carlo, based on the incoherent sum of Eq. flHU]) 
and Eq. (!90l) would possibly be less efficient than one based on (191]) . At the level of 
double-emission taken alone, this observation is of not much interest of course. 

Finally, let us mention that if the large limit is used, the interference between two 
distinct color contributions to (PTj) in the amplitude squared and averaged over spin degrees 
of freedom is suppressed by a factor A*"^. It is a direct consequence of the fact that 

E Tt{T^T'T^T') = Z^(l-l^ (93) 



a,b 

and 



J2 Tr(T''T^T''T'^) = ^ |^1 - — J 



(94) 



That completes the study of this kinematical configuration. We have explicitly trun- 
cated the complete amplitude to obtain an expression consisting of two incoherent dipole 
chains, which are easy to implement into a Monte Carlo program incorporated with exact 
Lorentz invariant phase space as in QED [1 , 2] for the solution based on the exact matrix 
element and full phase space coverage. 

Note that the truncated amplitude is given in the same kinematical formulation as the 
complete exact amplitude for two-gluon emission. It can be calculated at any point in phase 
space, and not only where the approximation is justified. That is why the re-installation of 
the exact amplitude into a calculation based on sufficiently refined dipole language can be 
accommodated for by a well-defined weight, namely the ratio of matrix elements squared. 
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We can conclude that this part of the exercise is finished as far as the analysis of two- 
gluon emission amplitudes is concerned. It can serve as a basis for studies of algorithms 
where the soft-gluon approximation is used as first step in event construction. The next 
step of our work must include other processes of double-parton emissions taken at tree 
level. Only then we can be sure that the method as used in the PHOTOS Monte Carlo 
can be of any interest for QCD. We want to note, that the picture of linked dipoles is used 
with great success in Ariadne parton shower [40] since quite some time now. 



5.1.1 Two soft gluons, general case 

Let us modify the previously discussed case now. We will demand that i/i ^ k^{oT A;^) 
but we will not ignore regions of phase space where k^k2 3> ki-k2 (this condition is to 
be taken again in the reaction frame). Under such conditions we may neglect all terms 
like {^il/ti)/{2p-ki). If the scalar product ki-k2 appears in a propagator denominator to- 
gether with larger terms, it can be dropped from that denominator. Thanks to the above 
approximation, we find that (l55ll . (l56l) . ( |571) and (1581) simplify again 



(1,2) 
mix 



p-ei k2-ei\ f ki-e2 q-e2 



p-ki k2-ki J \ki-k2 q-k 



<-2 



^ ^ p-k2 f p-ei k2-ei \fp-e2 fci-ea ', , - 

+ / 



p-ki+ p-k2\p-ki k2-kij\p-k2 ki-k2 
q-ki fq-ei k2-ei\ f q-e2 ki-e2 



q-ki + q-k2\q-ki k2-ki J \q-k2 ki-h 



_ J P-(^2 ki-e2\(k2-ei q-Ci 



mix 



p-k2 ki-k2j \k2-ki q-ki 
+ (96) 



+ / 



p-k2+ p-ki\p-k2 ki-k2j\p-ki k2-ki 
q-k2 f q-e2 ki-e2\ f q-Ci k2-ei 



q-k2 + q-ki\q-k2 ki-k2 J \q-ki k2-ki 



j[i,2] = \ ^( P-^i -y-h ^ q-k2 - q-kA f fca-ei fci-ea _ ea-ei , ^^^^ 



2 \p-k2+p-ki q-k2 + q-ki J \k2-ki k2-ki k2-ki^ 

Both /^jf'' and I^-^ are now build from scalar factors and /. The /jj^-^^-part, potentially a 
contribution to running coupling constant, remains this time. Recall that it is free of sin- 
gularities when the gluons become coUinear thanks to a cancellation between contributions 
of the quark and anti-quark lines. All other parts, and in particular all parts of /^ix ' could 
be dropped out as non-dominant. Let us study the terms surviving the approximation. 
After minor manipulations, we find that the approximated amplitude can be written as 
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the sum of two terms, the first of which is given by 



-v{p)Jfu{q) 



X 





p-k2 




_p-ki +p-k2 


\p-ki 





k2-ei\fp-e2 q-e2 



+ 



q-ki 



k2-ki 
'q-ei 



I rpbrpa 



p-ki 



q-ki + q-k2 \q-ki 
p-e2 ki-e2 



p-k2+p-ki 

+ 



p-k2 
q-k2 



ki-k2 
'q-e2 



p-k2 

p-ki 
p-ei 
p-ki 

p-e2 



q-k2 
q-e2 



q-k2 
q-ei 



q-k2 + q-ki \q-k2 p-k2 ) \q-ki k2-ki 



q-ki 
q-ei 



ki-e2 
ki-k2 



k2-ei 



+ 



p-ki - p-k2 ^ q-k2 - q-ki 



X 



p-ki +p-k2 
ki-e2 f p-ei 



+ 



q-ki 

q-ei 



krk2 \p-ki q-k 



q-k2 

+ 



p-e2 _^ q-e2 



p-k2 q-k2jk2-ki 



k2-ei 



p-e2 q-ei _ p^ g-ea 
p-k2 q-ki p-ki q-k2 



62-61 
k2-ki 



(98) 



The expression differs from the one discussed in the previous subsection. More terms, of 

the linked dipole type, survive. They are weighted by spin-and color-independent factors. 

Some remnants from and I^^^ which survive the approximation are combined with 

[1 21 

the surviving pieces from P^^^. Together, these parts form the contribution to the coef- 
ficient of [T°-,T\ for the contribut ion to the running coupling constant. However, such 
a manipulation may considered to be somewhat supeficial. Note that it concerns terms 
proportional to 

^p-ki+p-k2 q-ki + q-k2j ' 
which approaches zero when ki and k2 become coUinear. It is of course essential in this 
that we do not separate the expression into a p-and a g'-part. The second term is given by 



p-ki 



p-k2 _^ q-k2 



q-ki 



X 



rj-iarj-ib 



rpbrpa 



p-ei q-ei 



p-ki 
p-e2 _ 
p-k2 

p-ei 



p-ki+p-k2 q-ki + q-k2 
ki-e2 p-e2 



q-ki 

q-e2 

q-k2 
q-ei 



k2-ki 

k2-ei 
ki-k2 
p-e2 



p-ki q-ki J \p-k2 q-k2 



q-k2 
p-ei 
q-ki 
q-e2 



+ 



q-ei 


k2 


eA 


f p-e2 


q-ki 


ki 


k2) 


\p-k2 


q-e2 


ki 




f P-ei 


q-k2 


k2 


ki) 


\p-ki 



q-e2 
q-k2 

q-ki 



(100) 



It is again non-leading but can be interpreted as a correction to the dipole picture. 

Finally let us observe that it is straightforward to reproduce the amplitude of the 
previous subsection. If we would assume that either ki <^ /c2 or /c2 -C ^1, then all spin and 
color independent factors reduce either to zero or to one. 
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5.2 pt ordering or DGLAP 



In this subsection, we use the type of phase space organization (or in fact approximation) 
that is called Pt ordering. As in the previous case, it can be formulated in terms of 
properties of Lorentz invariant expressions, allowing for an easy identification of the parts 
of the amplitudes which may be simply dropped out. In the following, we will discuss the 
two kinematical cases of emissions into one, or two (opposite) hemispheres. We will devote 
a separate subsection to the discussion of the running coupling constant contribution. 



5.2.1 Dominant parts for emissions in one hemisphere 

We assume first that p-ki p-k2 or q-ki ^ q-k2, but we accept configurations where 
p-ki ~ ki-k2 or q-ki ~ ki-k2. In practice, such conditions means that k2, p and q are 
basically parallel to each other from the point of view of the direction of ki. The 'macro- 
scopic' size of ki-k2 implicitly removes the phase space regions where the gluons become 
coUinear. Such a case must be treated separately. Under our conditions, the T'^T^ part of 
the amplitude takes the form: 

j(i,2) ^ f _ k2-ei _ ^1^1 \ ^2^2 ^1-62 _ q-e2 
~ \p-ki k2-ki 2p-kJ^\2q-k2 krk2 q-k2 

q-h fq-ei k2-ei \ /g-e2 ki-e2 . ^^^^ 

q-ki — ki-k2 \q-ki k2-ki 2q-kiJ\q-k2 ki-k2 2q-k2 

Analogously, the T^T"' part reads: 

^(2,1) _ f _ ki-e2 _ jih \ J hh _^ ^2-ei _ g-ei 
\p-k2 ki-k2 2p-/c2/ \2q-kx k2-ki q-ki 
^ p-ki fp-e2 ki-e2 ^^2 Wp-ei k2-ei ■ 



mix 



p-ki — k2-ki \p-k2 ki-k2 2p-k2j\p-ki k2-ki 2p-ki 
It is only a matter of a short manipulations for the T'^T^ amplitude part to take the form: 

^(1,2) _ ki-k2 f p-ei _ k2-ei _ 4ih_\ J hh ^ ^1-62 _ g-e2 
q-ki - ki-k2\p-ki k2-ki 2p-kiJ \2q-k2 ki-k2 q-k2 



mix 



q-ki f q-ei k2-ei ^1^1 \ f p-ei k2-ei 



q-ki — ki-k2{ \Q-ki k2-ki 2q-kiJ \p-ki k2-ki 2p-ki 



and analogously, for the T^'T" part we find: 



^(2,1) _ ki-k2 f p-e2 ki-e2 i2h\ f[ , ^2-ei q-ei 



mix 



+ 



p-ki — k2-ki\p-k2 ki-k2 2p-k2j \2q-ki k2-ki q-k\ 
p-ki f p-e2 ki-e2 ^2^ 



p-ki — k2-ki \p-k2 ki-k2 2p-k2 

p-ei _ k2-ei _ ^ 

p-ki k2-ki 2p-kiJ^ ^\2q-ki k2-ki q-k 



.,P:e,_A^_^\ ^ /J^^A^_g^- ^^^^^ 
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The above expressions can easily be understood in the language of evolution kernels. 
The factors (q-ki)/ {q-ki — k2-ki) and (p-ki) /{p-ki — /c2'^i) can be understood as redef- 
inition of the spinor normalization, again similarly as explained in [33]. Note that the 
contents of the curly brackets in the last two formulas are free of terms proportional to 
{ki-e2)/{k2-ki) and {k2-ei)/{k2-ki). They represent simple dipoles spanned on the p,g-pair. 
In contrast with the previous section, correcting terms with I/ti and/or ^2 in numerator re- 
main. This is closely related to the necessity/possibility to introduce effective momenta 
flowing into J, once the language of PDFs is introduced. Note, that thanks to the pt 
ordering there are only two chains of dipoles left, similarly as it was in the BFKL-case. 

This time the coefficient (1631) for the running coupling constant seems to survive in its 
near complete form. In our limit and because of cancellations it develops an extra power 
of ki-k2 in the eikonal part, and a factor of \/ki-k2 in the terms proportional to ki^2- This 
will cancel the singularity of the virtual gluon once the amplitude is squared and a partial 
integration over the phase space is performed. Because of this, an interpretation in the 
language of the running coupling constant is prevented. It can be treated as a non-singular 
correction, and thus ignored. At this point, the freedom of choice for the arguments of the 
running coupling constant remains unrestricted. Note also, that remnants from the dipole 
parts did not migrate into the coupling constant expression, in contrast with the situation 
for BFKL in the previous section. This observation can be justified further. If we return 
to the soft gluon case of the previous section, and apply the pr ordering as well, then the 
running coupling constant contribution of Eq. (1981) will become suppressed. We will return 
to this point later in Section I5.2.3[ 

For now, we can conclude that our amplitude, with properties consistent with pt or- 
dering, is more compact than the exact one. It is gauge invariant and valid all over phase 
space. Such a simple, amplitude-level, expression can be used at the intermediate step for 
the definition of a parton shower algorithm, or to better understand the already existing 
ones. As in the previous BFKL case, an explicit form for the weight necessary to reinstall 
the exact distribution based on the two-gluon amplitude is available. It was mentioned al- 
ready earlier how to identify emission kernels already at the spin amplitude level. It seems 
that the explanations included in [33] can be applied here as well. For a further continu- 
ation of our discussion, the language of the truncated exact spin amplitudes is probably 
too restrictive: a particular choice for the direction of ki and/or k2 in the vicinity of p 
(or q) must be made, and generates a multitude of options. Furthermore, a definition of 
intermediate states for effective fermions valid all over the phase space would be necessary. 



5.2.2 Dominant parts for simultaneous emissions from p and q 

Now, let us assume that p-ki ^ p-k2 and q-k2 ^ q-ki. Under such conditions, k2 is ba- 
sically parallel to p, and ki to q. It is obvious that contribution PHI) dominates over all 
other contributions of the complete spin amplitude. If in addition p-ki ^ ^'■^2, then we 
can replace {k2-ei) / {k2-ki) with {p-ei) / {p-ki) and get 

j(2a)_/P-^2 ki-e2 hh \ J hh ^ P-ei /^q^x 
\p-k2 ki-k2 2p-k2) \2q-ki p-ki q-ki) 
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This expression is again of the dipole form. At the same time it is part of the expression dis- 
cussed previously. This means, that the replacement of {k2-ei) / {k2-ki) with (p-ei)/(p-A;i), 
which is potentially dangerous and valid only in this particular region of phase space, is not 
necessary, except for our proof here. Difficulties, as the ones to be discussed in Section 15.31 
can be avoided. 



5.2.3 Case when ki and k2 may become parallel 

So far, we have implicitly excluded the phase space region contributing to the running 
coupling constant from our discussion. The px ordering makes it unfavorable for the two 
gluons to become parallel one to another. Therefore, we have to look at regions of phase 
space where the virtuality of the gluon may approach zero separately. For that purpose, we 
will assume that the overall pt of the virtual gluon is small, but larger than ki-k2. We have 
found such an approach in studies presented in [41,42]. We will consider the configuration 
when q-ki 3> p-ki, q-k2 ^ ^-^2, p-ki ^ ki-k2 andp-fci ~ p-k2. Such choice represents the 
splitting the anti-quark with momentum p into a virtual fermion line entering J, and a 
single gluon of small virtuality and (moderately) small pt- Under such circumstances, 
dominant contributions from (1551) . ( |56l) . ( l57j) take the form: 

j(i,2) _ / p-ei _ k2-ei _ \ J ^^2 ki-e2 _ q-e2 \ 
"'^ ~ U-^i h-ki 2p-ki)\2q-k2 k^-k2 g-fcs/ 
_^ p-k2 f p-ei _ k2-ei _ / p-e2 _ ki-e2 

p-ki + p-k2 \p-ki k2-ki 2p-ki)\p-k2 ki-k2 

+ j—^ih— ( _ _ \ f _ ki-e2 

q-ki + q-k2 \q-ki k2-ki 2q-ki)\q-k2 ki-k2 




(106) 



2q-k2j 



7-(2,l) 
mix 



p-e2 ki-e2 ^2h \ J , ^2-61 q-ei 



p-k2 ki-k2 2p-k2 J \2q-ki k2-ki q-ki 
^ p-ki f p-e2 ki-e2 hh \f p-(^i h-ei ^1^1 ^^^^^ 



+ 



p-k2+ p-ki \p-k2 ki-k2 2p-k2j\p-ki k2-ki 2p-ki 
q-k2 f q-e2 ki-e2 |^2^2 A /g-ei k2-ei ^i^i 



q-k2 + q-ki \q-k2 ki-k2 2q-k2j\q-ki k2-ki 2q-k 



j[i,2] _ i p-ki-p-k2 ^ q-k2 - q-ki \ / fci-ea k2-ei _ 61-62 
~ 2'^\p-ki+p-k2 q-ki + q-k2j\ki-k2krk2 ki-k2 
1 1 [hhhh-hhhV^ 



Ap-ki+p-k2\ ki-k2 
1 » 1 / hiihh - hhhii 



4 q-ki + q-k2\ ki-k 



Contribution (155]) becomes non-leading as a whole, so we have dropped it out already 
now. Obviously our choice of remaining terms is different from the one of Section 15.2. 1[ 
Nonetheless one can observe that (11061) and (I107p are included in (11031) and (11041) . if for 
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the latter both possibihties p-ki ^ p-k2 and p-ki <^p-k2 are added together. Expression 
(11081) may contribute to the running couphng constant. The last line of (11081) is non- leading 
for our choice of kinematical conditions and should be dropped out. However, it is gauge 
invariant and rather compact, and it is necessary in case gluons would be collimated with 
q. Also, one should bear in mind that the first line of (I108P is less singular than it may 
seem, because of a partial cancellation of p and q dependent terms, which one may be 
tempted to treat separately. This is again, as in BFKL case, a consequence of the gauge 
cancellation, which is not exact as it concerns a virtual gluon. 

5.3 Angle ordering or CCFM style 

The CCFM case is definitely more difficult for our approach than the ones discussed so 
far. We have to use angular ordering to select the dominant terms. Such a choice is 
less straightforward for the Lorentz invariant representation. It is thus not possible to 
simply neglect some terms because they would be explicitly smaller than others. In fact, 
it seems that all terms of our expressions for the amplitude will need to be kept. Some 
kind of a language exploiting the concept of effective intermediate states is needed, and 
that is definitely a complication. Even if approximated amplitudes would be defined, they 
would differ from the complete ones quite substantially by the presence of these effective 
intermediate states. That would definitely be out of scope of the present paper. 

On the other hand, we want to remark that the interpretation of the result using 
some sort of the dipole language will persist in this case as well, as it is already visible 
in the exact amplitudes. There will be simply more terms to take into account, that is 
formulas (1431) . (I44p . fHS]) . (1491) . (I50p . The fermion normalization redefinition factors: 

{q-ki)/{q-ki + q-k2 — k2-ki) and {p-ki)/{p-ki + p-k2 — k2-ki) are more complicated than 
in the previously discussed cases. Also the terms contributing to the running of the coupling 
constant would require a more sophisticated treatment. 

6 Summary 

In this paper we have analyzed different forms of the exact tree level QCD spin amplitude 
for the process qq Jgg where J may represent any, color singlet, current. We have 
found quite well structured representations consisting of sums and products of compact 
gauge invariant parts. Some of those parts have a well defined meaning, as dominant 
contributions necessary to define evolution kernels. In contrast to previous studies for 
QED, the present expressions do not seem to be unique. This is understandable, since 
amplitudes and the structure of singularities are more complex in QCD. In particular, the 
discussion of singularities for processes where incoming quarks would be interchanged with 
outgoing gluons should be included in our considerations to constrain ambiguities. 

We have used two types of organizations of parts. The first one, discussed in Section [31 is 
rather directly obtained from Feynman rules, with only a limited set of simplifying algebraic 
manipulations. This form can be useful, because it manifests the relation between QCD 
and QED amplitudes. In Section |U our attention turned towards a form useful for QCD 
phenomenology. Our final expression consist of a sum of 10 gauge invariant parts, each just 
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one line long. Even at the level of the exact spin amplitude, one could attempt to provide 
a physical interpretation for these parts. We have provided such an attempt in Section H31 
The results look encouraging, and provide a separation into terms possible to interpret 
either as responsible for consecutive emissions of gluons or as contribution to the decay of 
a virtual gluon (running of the coupling constant), they however are superficial to a certain 
degree. In approximations, a migration of terms takes place between the gauge invariant 
parts as organized in the exact amplitude. Such a phenomenon should be expected and it 
is probably related to the choice of scales for the running coupling as accommodated by 
different schemes of parton shower models [43,44]. 

In Section [5l we discussed approximations for the amplitude, which consist of BFKL 
or DGLAP pictures, valid in the appropriate phase space regions. As a consequence, 
some parts of our expressions for the exact amplitude could be dropped and even shorter 
expressions appeared. They are straightforward to interpret. At the same time explicit 
expressions for the difference with the exact amplitude at any point of phase space exist. 
Such correction is straightforward to manipulate in any numerical application, for example 
in parton shower Monte Carlo applications, using positively defined weights. 

The discussion of simplification presented here is by far incomplete with respect to 
other discussions present in the literature. For example, we have not performed partial 
integration of the phase space variables nor have we introduced intermediate effective quark 
fields into the description. In particular, our discussion covered the CCFM picture only 
marginally, even though we do not neglect any aspects related to coherence. 

The main analytical results are the compact expressions for the amplitude given in 
Section [3] and Section 14.21 The discussion of approximations is fragmentary only, however, 
the localization of gauge invariant parts of the amplitude proved to be useful for that 
purpose. For a more complete discussion, we need to collect amplitudes for at least all 
processes of four external color fields accompanying the color- neutral current. 
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